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ABSTRACT 

Using advanced numerical schemes and grid-refinement we present 2D high-resolution 
models of solar granulation with particular emphasis to downflowing plumes. In the 
high-resolution portion of our simulation, a box measuring 1.97 x 2.58 Mm 2 (vertical 
x horizontal), the grid size is 1.82 x 2.84 km 2 . Calculations at the resolution usually 
applied in this type of simulations amount to only a few horizontal gridpoints for a 
downflowing plume. Due to the increased number of gridpoints in our high resolution 
domain the simulations show the development of vigorous secondary instabilities of 
both the plume's head and stem. Below a depth of about 1 Mm the plume produces 
patches of low density, temperature, pressure and high vorticity which may last for all 
of our simulation time, ~ 10 minutes, and probably considerably longer; they may be 
ascribed to the 2D nature of the present calculations. Centrifugal forces acting in these 
patches counteract the strong inward pressure. Probably most importantly, the plume's 
instabilities give rise to acoustic pulses created predominantly down to ~ 1.5 Mm. The 
pulses proceed laterally as well as upwards and are ubiquitous. Ultimately most of them 
emerge into the photosphere. A considerable part of the photospheric 'turbulence' in 
these models is due to those pulses rather than to some sort of eddies. - The upflows in 
granules are smooth where they reach the photosphere from below even in the present 
calculations; however, the pulses may enter in the photosphere also in granular upflows. 
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1 INTRODUCTION 

Considerable activity is currently going on aimed at observ- 
ing the solar atmosphere in high resolution. This is true for 
ground based observations (e.g. the Swedish vacuum tele- 
scope in Teneriffe), balloon borne instrumentation as in the 



Sunrise project, see Schiissler & The Sunrise Team (20031 



and satellite observations (e.g. the recently started HINODE 
satellite). The observations performed or planned with such 
facilities will contribute to a better understanding of hydro- 
dynamic and magnetohydrodynamic phenomena near the 
Sun's surface. It is mandatory that these observational ef- 
forts be paralleled by simulations of truly high resolution, 
having in mind ultimately the same goal. 

We want to make a step in this direction. Current sim- 
ulations of solar granulation and the atmosphere above it 
typically seem to apply a (horizontal) grid spacing down to 
about 25 km, e.g. |Stein fc Nordlund| ( |T^8j ) or |Steffen| ( [2007l |. 
Naturally, 3D simulations prevail nowadays. Since with de- 
creasing mesh size the computational load becomes quite 
soon prohibitive it is tempting to refine the computational 
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grid in a localized region of interesting processes only, using 
a normal grid spacing elsewhere. The technique of grid re- 
finement is quite customary in other areas of astrophysical 
simulations, for example in supernova research as in |Plewa| 
|fc Mulle7] ( |2001| l, but does not seem to have previously been 
applied in the solar context. 

When it matters to achieve high resolution the numeri- 
cal method should also be chosen with that goal in mind. We 
have adopted high-resolution schemes of the essentially non 
oscillatory (ENO) type. An overview over ENO methodol- 
ogy can be found in Harten et al. (1997 1. With high accu- 



racy these methods can handle smooth flow, steep gradients 
and discontinuities at the same time. They can also be run 
largely without resorting to artificial diffusivities. 

Our specific approach has its pros and cons. The most 
basic contra is the restriction of our simulations to two spa- 
tial dimensions. While this can and must be considered 
a drawback, there are, on the other hand, major advan- 
tages coming along with this decision. We can, in partic- 
ular, easily attain a spatial resolution in the refined area 
which is out of reach in 3D; the present calculations adopt 
a grid size of 1.78 x 2.84 km 2 (vertical x horizontal) in the 
high resolution area. Compare this grid spacing with the 
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presently typical values mentioned above and the resolu- 
tion of the solar optical telescope at the HINODE satellite: 
0.25 arcseconds ~ 175 km. (To be fair, we have to keep in 
mind in such a comparison that a single computational grid 
point cannot represent a physical feature.) The 2D models 
can already give hints towards phenomena which might oc- 
cur in real solar granulation and can suggest ways to look 
at observations. They can aid in more meaningful decisions 
what precisely should be investigated in 3D simulations. 
Note in addition that, even if 2D and 3D calculations nat- 
urally do not yield identical results, they still exhibit quite 



some similarity, see the discussion in Asplund et al. ( 2000 I 



As a consequence, even today 2D investigations are con- 
ducted with specific goals in mind, e.g. in |Wedemeyer-B ohm 
(2005l. Of course, the ultimate investigations of this 



et al. 



sort will have to be done in 3D; it is conceivable, however, 
that some of the findings are such that they may be en- 
countered in 3D as well, in particular the results of section 
3.2. 

Naturally, such calculations should in the long run ad- 
ditionally include magnetic fields. The present calculations 
consider the purely (radiation-) hydrodynamic case. They 
thus apply predominantly to the quiet Sun. It can be con- 
sidered useful to perform both HD and MHD simulations in 
order to assess which phenomena must basically be ascribed 
to hydrodynamics and which ones to magnetohydrodynam- 
ics proper. 

In this paper, we direct our attention primarily to the 
downflowing plumes and the processes they excite. After all, 
these plumes are the most vigorous entities in solar granu- 
lation, have been encountered in any numerical model and 
are obviously a major ingredient to the overall dynamics. 
In addition, they are expected to be crucial in exciting so- 



lar p-mode oscillations, see Rast ( 1999 I and the references 



there. In many simulations addressing either solar granula- 
tion or more idealized models of compressible convection, 
these plumes penetrate down several pressure scale- heights. 
Deeper down, by possibly joining with other plumes in a 
specific manner, they give rise to mesogranulation, |Stein fc| 
Nordlund (1998 l, and perhaps even supergranulation. 



2 NUMERICAL METHODS AND PHYSICAL 
MODEL 

2.1 Numerical Methods 

We have developed the package ANTARES (A Numerical 
Tool for Astrophysical RE Search) for simulations of solar 
granulation and other astrophysical flows. A more detailed 



description of ANTARES will be given elsewhere (Muthsam 
et al.|2007 i 



In the mode which is appropriate here ANTARES al- 
lows the numerical integration of the continuity equation 
and the equations for momentum and energy balance in 2D 
(also in ID or 3D) in time. The equations are solved in con- 
servative form. The numerical scheme applied for the present 
calculations is the 5^ order accurate weighted ENO scheme 



of Liu et al. (19941. In fact, we have various, quite different 



ENO schemes implemented and checked against standard 
examples and against each other; [Muthsam et al.| ( |2007| |. - 
Time integration is achieved using specific second or third 



order Runge-Kutta schemes that preserve the properties of 
ENO schemes. 

We solve the equation of radiative transfer by the 
method of short characteristics, 



Mihalas et al. 



(1978) 



This method makes ample use of interpolations which are 
done to at least second order of accuracy. At high optical 
depth we switch to the diffusion approximation, as appro- 
priate. Values for opacity coefficients and the equation of 



state are taken from the OPAL tables, Rogers et al. (1996 l 
augmented by line opacities based on the 



Kurucz 



(1992) ta- 



bles; the data as we need them are generated by an extension 



of the Piskunov & Kupka (2001 1 code. Wavelength variation 



of the opacities is accounted for by the method of opacity 



binning according to Nordlund ( 1982 I, Ludwig et al. ( 1994 I 



This amounts to group into the same bin frequencies v for 
which optical depth t,, = 1 is reached at similar geometrical 
depths. 

Fixed, slip-free plates provide the boundary conditions 
at the top and at the bottom. At the bottom, the amount of 
energy flux appropriate for the solar case enters the domain 
via radiative diffusion. 

The fixed boundary conditions at the top lead, in some 
instances, to a reflection of pulses. Still, since they are ap- 
plied high in the atmosphere at low densities, these artificial 
effects are felt only in the higher atmosphere and not in the 
bulk of our domain. No major unwanted influence is exerted 
by the lower boundary condition. 

2.2 Physical and numerical parameters 

The basic computational domain without grid refinement 
encomprises a size of 2.78 x 11.19 Mm 2 (depthxwidth) at a 
resolution of 382 x 985 equidistant grid points. The dimen- 
sions of a numerical cell consequently are 7.28 x 11.36 km 2 . 
The density at the upper boundary of the whole domain is 
~ 2 x 10~ 9 g cm -1 . This corresponds to a location with a 
small optical depth in the continuum (~ 10 -4 ) according 
to any standard solar model atmosphere. The basic model 
was relaxed to a flux constancy of ~ 1% and then the high- 
resolution part was plugged in. 

The refined region measures 1.97 x 2.58 Mm 2 with 
1081 x 909 gridpoints, yielding a size of 1.82 x 2.84 km 2 for 
the numerical cell. The high resolution domain starts also 
quite high in the atmosphere and spans ~ 7 pressure scale 
heights downwards. 

Radiative transfer makes use of 12 bins. We switch 
smoothly from full radiative transfer to the diffusion ap- 
proximation at an optical depth of tr 03S ~ 10000. 

Regarding the diffusion of momentum we have per- 
formed two sets of calculations. In model 2 a very small diffu- 
sion term is included in the momentum equation (only) such 
that for the Prandtl number Pr we have formally Pr < 0.01 
everywhere. In model 1 we have added diffusion terms in the 
momentum and energy equations according to the artificial 
Vogler et al. | ( |2003) , who again re- 



diffusion description as in 



fer to Stein & Nordlund (1998 l. This procedure amounts to 



adding a hyperviscous term and a term aimed at broadening 
of shocks to the equations. This leads to an effective Prandtl 
number still quite small in the upper regions of the simulated 
domain; however, for model 1 the effective Prandtl number 
gets O(l) in the lower two thirds of our computational box, 
in particular at places of strong velocity gradients. - On a 



High-resolution models of solar granulation 3 



quantitative level, these two models differ to some extent. 
They agree, however, on the more qualitative level which is 
relevant for the present paper. Since model 2 has not been 
run quite as long as model 1 the presentation given here 
focusses on model 1. 



3 DISCUSSION OF THE SIMULATIONS 
3.1 The plumes 

3.1.1 The plume head 

The generation of new solar granules often commences with 
the appearance of a more or less circular spot of low tem- 
perature and high density forming through cooling effects 
in the middle of a preexisting granule. Due to the rapid de- 
crease of opacity with decreasing temperature this cooling 
constitutes a runaway process and initiates a downflowing 
plume of dense, cool material. In the real Sun, i.e. in 3D, 
the plume later on frequently extends sidewards along lanes 
and causes in this way a disruption of the original granule. 
The plume tends to be surrounded by more gentle upflows. 

Considering now the evolution of a plume according to 
our calculations, we start with a discussion of the plume's 
head. In Fig. 1, three main plumes are visible. Somewhat to 
the left of the center there is a nascent plume. The nascent 
plume's boundary and, by the way, also the lower boundary 
of the photosphere above granular upflow are exceedingly 
smooth and sharp even at this resolution, judging from spe- 
cific entropy, which is plotted in the figure, and also from 
other relevant quantities. This finding is in agreement with 
3D results obtained previously at common resolution as dis- 
cussed above (Stein & Nordlund 19981.- The next plume 



to the left is in an advanced stage of stem instability (see 
the next subsection). Still farther left again is a plume in 
an earlier stage of development, exhibiting head instability. 
This instability comes about in the following way. Very soon 
after plume formation there is a tendency for the head to 
be faster than the subsequent material (and to broaden at 
the same time). Sidewards of the following tail low-density 
regions may easily develop. Their material flows in part to- 
wards the rear side of the head and gets entrained by the 
head's sideward lobes. In the course of time these rapidly 
spinning regions of low density detach from the plume proper 
and remain in appreciable depth (> 1 Mm) for quite some 
time, about 10 minutes as much as we can tell and possi- 
bly considerably longer, although the contrasts to the sur- 
rounding material are really large, since its density may be 
only half of the outside value. These patches of high vortic- 
ity and low density resist pressure via Coriolis forces. They 
show no tendency to ascend and to appear in the photo- 
sphere within a granular upflow. For further discussion of 
these vortex patches see subsection 3.1.2. 

Instabilities of downward plunging plumes in compress- 
ible, layered media have been analyzed by[Rast (1998). The 
settings in that paper were idealized, the initial structure 
of the layer being related to polytropes. In particular, |Rast| 
( 1998 ) investigated the case of low and high Prandtl number 
Pr. For Pr small (Pr ^ 0.1) the head is a small, essentially 
circular, rotating structure which eventually detaches from 
the tail, with a subsequent similar structure appearing again 
and detaching from what is now the plume's head. On the 



other hand, the head's structure in the high-Prandtl-number 
case is much more similar to what we observe. Since we have 
a low Prandtl number case (using only the realistic radiative 
conductivity and, in some calculations, no artificial diffusiv- 
ity) the question about the cause of that discrepancy arises. 
For a resolution of that issue consider the following. 



Firstly, as judging from figure 3 of Rast (1998 l, the hori 



zontal extent of the initial disturbance assumed there seems 
to be quite small compared to the pressure scale height. In 
our case, the cool spot which develops initially in the orig- 
inal granule is quite broad, of the order of a pressure scale 
height; in addition, its horizontal temperature profile does 
not resemble a peaked Gaussian but is trough-like instead. 
Furthermore, we have a general velocity field present, which 
will, via shearing forces, efficiently act on this broad initial 
disturbance, facilitating the development of further insta- 
bilities. We conclude that for these reasons we do not nec- 
essarily have to expect agreement with the results of |Rast| 
(p98|. 



3.1.2 The plume stem 

Across the plume stem and in particular at the boundary 
to the surrounding material strong velocity gradients are 
present. As a consequence, as soon as the plume stem has 
attained some length a strong Kelvin-Helmholtz instability 
sets in. If surrounding disturbances are not overly large, the 
plume has a chance to initially develop a sinusoidal, lateral 
disturbance. Very soon, a von Karman vortex street appears 
which is well known to be a typical 2D instability. These 
vortices (as well as similar vortices which have been formed 
in the case of more severely disturbed plumes) are rapidly 
spinning, nearly circular patches of low density gas. As the 
plume disrupts due to increasing instability the vortices get 
separate entities, quite similar to the vortex patches which 
have been created out of the head's instabilities. 

Such vortex patches are already visible in our large, un- 
refined domain. Once the motions in our refined subdomain 
have become sufficiently turbulent, several of them are seen 
there at each instant of time. The truly persistent vortex 
patches are located in the lower half of our domain. 

In the context of idealized microphysics, similar phe- 
nomena have been observed by |Porter fc Woodward| ( |1 994| > 
in their simulation of 2D convection. These authors assume, 
in particular, a gamma-law gas and a constant heat con- 
duction coefficient. We agree in the solar context also on 
their observation, according to which density variation from 
the boundary to the center of a vortex patch is much larger 
than temperature variation. In a specific case we have found, 
for example, a density contrast of about a factor 2 between 
patch center and the ambient medium, whereas temperature 
changes by no more than ~ 20%. 

It should be noted, however, that those vortex patches 
are likely to be specific for the 2D case. They are not known 
in the 3D case. In particular, they do not seem to show up in 
our 3D high resolution models at least as far as these have 
been evolved up to now (i.e., not yet for a long physical 
time). In our 3D calculations we apply a grid refinement 
strategy similar to that one of the present paper and achieve, 
in the high resolution subdomain, a grid-size of 7.59 x 9.80 x 
9.80 km 3 (7.59 km = vertical); see 



Muthsam et al. 



(2007) 



As a consequence of all those instabilities, the undis- 
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turbed plume proper does not reach very deep (i.e., it usually 
ends, as far as we have observed, within our high-resolution 
subdomain). Still, plumes may nevertheless be embedded 
in broader, more gentle downflows, which reach essentially 
down to the lower boundary of our whole computational 
domain and would extend presumably even deeper if our 
domain would only extend farther down. Therefore, the in- 
creased instabilities of the plumes which we see in high 
resolution do not contradict the findings by Stein & Nord- 
lund ( 1998 I that downflows merge at large depth, forming 



the mesogranulation and possibly even the supergranulation 
network. 

As long as the plume's stem fairly maintains its integrity 
the velocities of the plume's core as measured in the fixed 
coordinate system tend to exceed the local velocity of sound. 

3.2 Acoustic pulses and atmospheric turbulence 

The downflowing plumes generate acoustic pulses. To a con- 
siderable part, pulse generation is associated with the insta- 
bilities of the plume's head and stem. For example, a stem 
may experience a strong lateral displacement, say to the 
right. Then the rapidly downflowing material can happen to 
collide with material belonging to a gentle upstream to the 
right of the stem thus causing a density, temperature and 
pressure enhancement. The location of disturbance genera- 
tion will, in general, move during the course of time. As the 
disturbance signals are continuously created at that moving 
point and propagate with sound speed they appear initially 
essentially as a straight or curved line, depending on sound 
speed variations and possibly also variations in the velocity 
of the source or the ambient medium. 

The initial size (length of a front) of a pulse is largely 
set by the source's lifetime. As a consequence, pulses gener- 
ated in the complex inflow field near the photosphere (where 
the pulses, in addition, easily escape upwards) may be local- 
ized; also, in the late stage of plume development where the 
integrity of the plume has been lost small or not so regular 
pulses may be created. As is evident from the generation 
mechanism described above these pulses are formed where 
flow fields of the necessary characteristics are present, i.e. 
mainly from somewhat below the photosphere down to a 
depth of ~ 1.5 Mm. Below about that depth the plumes 
are, in general, no longer sufficiently vigorous in order to 
generate pulses. 

Following the pulse front is a region of relatively high 
density and temperature, at least near its origin. Density and 
temperature are next to discontinuous across the pulse. The 
pulses can move against the local flow field and are therefore 
not advection phenomena but akin to sound waves. 

Due to strengthening mechanisms (see below) some 
pulses can travel over surprisigly long distances, some of 
them at least a few granule widths. - Pulses are ubiquitous. 
At each instance of time there is a considerable number of 
them in our high-resolution computational domain, in par- 
ticular also in the atmosphere, see Fig. 2, once the solution 
has only been given time to evolve away from the more lam- 
inar initial state taken from the coarse grid. Through inter- 
action with the complex inflow field of a granule the plumes 
may get enhanced just there so that chances are good to see 
very strong pulses (quite possibly more than one at a time) 
moving from below into the photosphere above a downflow, 



contributing to a very considerable part to turbulence (gra- 
dients of velocity etc.) there. They run, however, into the 
photosphere also in otherwise smooth granular upflows. 

The pulses, travelling with sound speed, exhibit in part 
linear, in part nonlinear character. Two oppositely moving 
pulses running into each other penetrate and continue their 
course more or less undisturbed, in the sense of wave super- 
position, which can also occur for wave packets. The same 
may also be true to some degree if a pulse traverses the more 
regular, not chaotic, part of a plume, although we believe to 
have seen plume strengthening as a consequence. 

Predominantly when a pulse hits a part of a plume with 
a complicated flow field, for instance in the lower part of 
a developed plume or in the atmospheric inflow region at 
the top, and interacts with the complex pattern of motions 
prevailing there, the interaction is nonlinear. In addition to 
the obvious distortion due to the variation in sound speed 
across a plume together with the effects of the plume's shape 
and flow field, it may trigger the generation of subsequent 
pulses or may be strengthened. 

Because of strengthening it may cross several downflow- 
ing plumes and maintain its identity over several granule 
widths and traverse, e.g., all the width of our high resolu- 
tion domain, > 2 Mm. Indeed, the dominant pulse in Fig. 2 
(the semicircle in the left upper part of the figure) is an 
example for the nonlocality of pulses. It has entered the 
high-resolution domain from the right as a quite feeble dis- 
turbance and gained in strength through repeated interac- 
tion with the complicated flow field of the plumes it crossed. 
Other pulses may be more local, among others because they 
have been directed mainly upwards right from the beginning. 

On the other hand, pulses do not seem to get easily 
absorbed when crossing plumes. We have seen, however, a 
reduction of the number of pulses procedding in the follow- 
ing way. A leading pulse leaves behind it a region of higher 
sound speed. When a trailing pulse is about to overtake the 
leading one these two pulses merge into one. - The pulses 
are noticeably refracted when entering the photosphere from 
below at some angle because of the sharp drop in sound ve- 
locity. 

The pulses and their generation mechanism obviously 
are not identical to the acoustic events discussed in lSkartlienl 



et al. (2000l. They found pulses which are generated when 



a small granule collapses and a larger downflow is initiated 
at the site of the collapse. Given the relatively rare starting 
configuration one may presume that these pulses are consid- 
erably less frequent than the acoustic pulses discussed here. 
- According to our evidence the collision of plumes, which 
seems occasionally to have been held responsible for the gen- 
eration of pulses, cannot be considered a major source of 
pulses at least from the viewpoint of frequency of events. 

These phenomena (secondary instabilities of the 
plumes, pulse generation) are most clearly seen in the high- 
resolution part of our calculation. It is remarkable that, how- 
ever, they are also distinctly present in the low resolution 
portion. We take this as an indication of the quality of the 
basic numerical scheme. 

It will be of interest to look into the material in a more 
quantitative fashion in order to figure out what these insta- 
bilities respectively pulses imply for the transport properties 
of convection and for questions of heating of the chromo- 
sphere in the quiet sun. Similarly, the possible role of the 
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acoustic pulses in the excitation of p-modes needs to be 
investigated. Of course, due to the main limitation of the 
present models, namely that they are 2D, it will be neces- 
sary to perform analogous simulations for the full 3D case. 
As mentioned previously, 3D calculations are presently being 
pursued. They apply a grid spacing of 7.59 x 9.80 x 9.80 km 
in the high resolution region, thus similar to the low reso- 
lution portion of the present 2D simulations, and will be 
reported in a forthcoming paper. 



ACKNOWLEDGMENTS 

This research was supported by the Austrian Science Foun- 
dation, project AP17024. We thank the referee for comments 
which helped to clarify the presentation. 



REFERENCES 

Asplund M., Ludwig H.G., Nordlund A., Stein, R.F., 2000, 
A&A 359, 669 

Harten, A., Enquist B., Osher S., Chakravarthy S.R., 1997, 
JCP 131, 3 

Kurucz, R.L., 1992, Opacities for Stellar Atmospheres: 
[+0.0], [+0.5], [+1.0]. Kurucz CD-ROM No. 2, Cambridge, 
Mass.: Smithsonian Astrophysical Observatory 
Liu X., Osher S., Chan T., 1994, JCP 115, 200 
Ludwig H.G., Jordan S., Steffen M., 1994, A&A 284, 105 
Mihalas D., Auer L.H., Mihalas, B.W., 1978, ApJ 220, 1001 
Muthsam H.J., Low-Baselli B., Obertscheider Chr., Langer 

M., Lenz P., Kupka F., 2007, in preparation 
Nordlund A., 1982, A&A 107, 1 
Piskunov, N., Kupka, F., 2001, ApJ 547, 1040 
Plewa T., Miiller E., 2001, Computer Phys. Comm. 138, 
101 

Porter D.H., Woodward P.R., 1994, ApJS, 93, 309 
Rogers F.J., Swenson F.J., Iglesias C.A., 1996, ApJ 456, 
902 

Rast M.P., 1998, Journal Fluid Mech. 369, 125 
Rast M.P., 1999, ApJ 524, 462 

Schiissler M., Sunrise Team, 2003, Proc. of the 16th 
ESA Symposium on European Rocket and Balloon Pro- 
grammes and Related Research (edited by B. Warmbein), 
ESA SP-530, pp. 279-283, ESA Publ. Div., Noordwijk 
Skartlien R., Stein R.F., Nordlund A., 2000, ApJ, 541, 468 
Stein R.F., Nordlund A., 1998, ApJ, 192, 91 
Stein R.F., Nordlund A., 2000, Solar Physics, 499, 914 
Steffen, M., 2007, p. 36 in 'Convection in Astrophys.', IAU 

Symp. 239, F. Kupka, I.W. Roxburgh, K.L.Chan (eds.) 
Vogler A., Shelyag S., Schiissler M., Cattaneo F., Emonet 
T., Linde T., 2003, p. 157 in 'Modelling of Stellar Atmo- 
spheres', IAU-Symp. 210, N. E. Piskunov, W. W. Weiss 
and D. F. Gray (eds.) 
Wedemeyer-Bohm S., Kamp I., Bruls J., Freytag B., 2005, 
A&A 438, 1043-1057 



H. Muthsam et al. 




Figure 1. Specific entropy (red = high) in the high-resolution domain. Note the sharp and smooth 
boundary of the nascent plume (left of the middle); the instabilities and entrainment of surface 
material (low entropy) in the more developed plume (farther left); and the head instabilites of the 
plume quite to the left. Compare with the smooth structures in the entropy plot in[Stein & No rdlund] 
(2000), fig. 11 there. In vorticity indications of turbulence are, however, seen there, I.e., fig. 13. 




Figure 2. Snapshot of the high resolution region with the plumes in a more developed state. Color: 
log (p/p) (p = density, p = horizontal average of density; red = high density. Isolines of log (p/p) 
(black), where p = pressure, crowd near acoustic pulses and shocks. The left half of the large pulse 
(semicircle) moves to the left, the other pulses crossing it there move to the right. A number of 
circular patches of low density and pressure and high vorticity can be discerned in the lower half of 
the figure. 



